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Abstract 



Single electron tunneling into small superconducting islands is sensitive to 
the gap energy of the excitations created in the process and, hence, depends 
on the electron number parity. At low temperature the properties of the sys- 
tem are 2e-periodic in the applied gate voltage, turning e-periodic at higher 
temperature. We evaluate the tunneling rates and determine the probabili- 
ties for the even and odd state as well as the cross-over condition from the 
balance between different processes. Our analysis includes excitations in the 
leads. They are essential if the leads are superconducting. The influence of 
parity effects on single electron tunneling, Cooper pair tunneling, and the 
Andreev reflection in superconducting transistors yields rich structures in the 
I-V characteristic, which explains recent experimental findings. 

PACS numbers: 73.40.Gk,74.50.+r,73.40.Rw 
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Recent experiments showed electron number parity effects in small superconducting is- 
lands |l|-|3|]. In an electron box, if the energy gap is smaller than the charging energy, the 
voltage shows a characteristic long-short cycle, with a 2e-periodic dependence on an applied 
gate voltage. This has been interpreted as the even-odd asymmetry predicted by Averin and 
Nazarov |4]]. At higher temperatures both parts of the cycle become more equal in length. 
Eventually above a cross-over temperature T cr the e-periodic behaviour typical for normal 
metal electron boxes is recovered. 

The even-odd asymmetry arises since single electron tunneling from an initial state with 
an even number of electrons on the superconducting island leads to a state with one extra 
electron - the odd one - in an excited state. In the small island, where charging effects prevent 
further tunneling, the odd electron does not find another excitation for recombination. As 
a result the energy of the final state lies above that of the equivalent normal system by 
the gap energy. Only at larger gate voltages another electron can enter the island, and 
the system relaxes to the ground state. The observed cross-over at higher temperatures to 
an e-periodic behavior was explained in Refs. |l]-[3| by free energy arguments for the island. 
This approach cannot be extended to describe nonequilibrium excitations in the leads. They 
are essential, however, if the leads are superconducting and are naturally included in the 
description presented below. 

We evaluate the transition rates due to electron tunneling between a small superconduct- 
ing island and normal or superconducting lead electrodes. The rates depend on the energy 
difference before and after the tunneling, which includes the gap energy Aj of the quasipar- 
ticle excitations in the island, and A; if the lead is superconducting. At low temperatures 
processes which cost energy are suppressed. This makes the tunneling of the one excited 
electron, the odd one, an important process. The rate of this process is smaller by a factor 
1/iVg than the total rate of the other N$ electrons in the condensate. On the other hand, 
the excitation energy in the island is regained. The balance between the rates determines 
the relative probabilities of even and odd number states. 

We recover the results of Ref. || in the case of a normal-superconducting (NS) electron 
box. We further provide the extension to the superconducting-superconducting (SS) electron 
box and superconducting transistors. The single electron tunneling current through SSS 
and NSN transistors shows at low temperature a 2e-periodic plateau structure (for suitable 
parameters a double plateau). In NSN transistors with typical parameters this current is 
several orders of magnitude higher than the cotunneling current studied in Ref. |3[]. At 
T > both the even and odd states are occupied. Then the Cooper pair tunneling in SSS 
transistors and the 2e-tunneling via Andreev reflection in NSN transistors also connect two 
even or two odd states. This results in two e-shifted, 2e-periodic sets of current peaks, 
recently detected in Refs. |]J5|||. 

At higher temperature the odd and the even state are equally populated leading to an e- 
periodic behaviour. The cross-over temperature T cr is of the order of the gap energy Aj + A/, 
divided by the logarithm of an effective number of electrons. The latter is large, resulting 
in a substantial reduction. 

We first consider an electron box , a closed circuit consisting of a small superconducting 
island connected via a tunnel junction with capacitance Cj to a lead electrode and via a 
capacitance C g to a voltage source U g . The charging energy of the system depends on 
Q g = CgUg and the number n of charges on the island 
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E ch (n,Q g ) - — . (1) 

The total capacitance is C = Cj + C g . The normal state conductance of the junction 
1/Rt = 47re 2 iVj(0)V;iV;(0)V;|T| 2 //i can be expressed by tunneling matrix elements. The 
island has the normal density of states iVj(O) (per unit volume) and volume Vi (similar for 
the lead). After a tunneling process, starting from an even state, one electron — the odd 
one — on the island is in an excited state. Therefore, the energy of the final state differs 
from the initial one not only by the charging energy but also by the energy of the excitations 
e « > Aj (and €\ > A; if the lead electrode is superconducting). 

Since the behaviour of the system is 2e-periodic in Q g it is sufficient to consider < 
Q g < 2e. We further assume that for < Q g < e the superconducting island is the even 
state. The transition rates due to single electron tunneling between the states with n = 
and n = 1 electron charges (up or down) on the island are []7||| 

This is the total rate due to the tunneling of one of the many electrons in the junction. It 
depends on the difference in the charging energy 5Ef h {Q g ) = ±[E ch (l,Q g ) — E ch (Q,Q g )]. 
The minimum energy difference is SE^ = ±5E c h(Q g ) + A» + A;. The energy of the exci- 
tations created in the tunneling process enter through the density of states into the well 
known classical quasiparticle tunneling characteristic it(eV) ||. The asymptotic form of 
the tunneling rates for T <ti A, in the SN junctions is 

r ( ^) = ("^ ex P(-f) for ^ >0 (3) 

I 7 ^ y /5E(5E-2A i ) for 5E < 0, 
For SS junctions we have for \8E\ > T near the onset 

' ' " I ^ fa SE < <4) 

The expressions (@-§) describes single electron tunneling from a state with even electron 
number to one with an odd number (more precisely from an initial state, where no electron 
is in an excited state). On the other hand, if we consider the tunneling from an odd electron 
number state (where one electron is in an excited state) we have to take into account two 
types of transitions: 

(i) One of the many electrons in the condensate of the island can tunnel back. The rate is also 
given by eq. (0). The minimum energy difference is 5E~(Q g ) = —SE ch (Q g ) + Aj + A;. Also 
in this transition we create two excitations, one on the island and one in the electrode. Hence 
5E~(Q g ) is positive in the interesting range of Q g , and at low temperature the transition rate 
is exponentially small. After such a process the two excitations in the island can recombine 
fast and return the island to the ground state. The extra excitation in the electrode will 
diffuse rapidly further into the leads and recombine there. Nevertheless, energy was needed 
to create the two excitations, which enters into the transition rate. 
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(ii) In the odd state there exists a second channel: the odd electron can tunnel back. Its 
initial energy is distributed thermally with the constraint J °° de^i^ei) f (ei) = 1. The 
tunneling rate of the odd electron into any one of the states in the lead is 

/■OO /"OO 

7 (Q a ) = 27r|T| 2 jV i (0)V i / de t / de^e^jV^) 

JO J —oo 

](e i )[l-f{e l )\5{e l -e i -5E ch ). (5) 



Here, J\fi/i(e) = 6(|e| — A.j/i)|e|/^/e 2 — A?^ are the reduced BCS densities of states. The 
rate is substantially different when the odd electron can tunnel from the lowest energy state 
6i = A, and when it has to be in a higher excited state (with exponentially low probability). 
Accordingly, 

J^dMT^S^ ior5E ch -A l + A i>T 

l(Q 9 )-\ V(A,/AQ Al _ SEch , A.-A,-^^ fnr A._A._J!R.NT (6) 



2e^ t JV i (0)V i ^( Al -SE ch r-A: 



: exp 



The rate 7 is reduced by a factor l/iVj(0)T^ compared to the transition rate T On 
the other hand, 7 may still be larger than T since it depends on the difference of the gaps 
5E = Ai — Ai — 5E c h, whereas the sum of the gaps enters into T. The three transitions 
between different states are visualized in the energy scheme in Fig. 1. Three states appear 
to play a role, but since the state with 2 excitations relaxes quickly to the ground state by 
quasiparticle recombination, effectively only one even state and the odd state need to be 
considered. However, the transition rate T oe from the odd to the even state is the sum of 
the two relevant rates 

T e °(Q g ) = T[5E + ] , r oe (Q 5 ) = T[8E-] + 7 • (7) 

The master equation for the occupation probabilities of these states W e (Q g ) and W (Q g ) 
is dW e (Q g )/dt = -T eo (Q g )W e (Q g ) + T° e (Q g )W (Q g ) with W e (Q g ) + W (Q g ) = 1. The 
equilibrium solutions are 

W e{o) (Q g )=T°< e °\Q g )/Tx(Q g ), (8) 

where r s (Q g ) = r oe (Q g ) + r eo (Q g ). For T oe > T eo we have W e (Q g ) w 1, i.e. the system 
occupies the even state, while for r eo 3> T oe the island is in the odd state. The cross-over 
occurs at W e ~ W a , i.e. 

r oe (Q cr )^r eo (Q cr ). (9) 

At low temperatures the odd electron tunneling rate 7 is large compared to T[5E~ (Q g )]. 
In this regime we find the even-odd asymmetric, 2e-periodic behaviour. At high temperature 
7 can be neglected and the cross-over between even and odd states occurs at the symmetry 
points Q cr = e/2, 3e/2, leading to an e-periodic behaviour in Q g . The cross-over tem- 
perature T cr between both regimes can be found from the condition 7(e/2) w T[5E + (e/2)]). 
Up to numerical coefficients inside the logarithm we find for T< Aj 
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T "-\£t s for A, < A, < 10) 

where N(T) = N (T) for A, < T < Aj and JV(T) = iV 5 otherwise. The numbers N (T) = 
Ni (0)Vi\ / 2n A iT and iVg = 7riVi(0)VjAj are the number of states available for quasiparticles 
near the gap and the effective number of superconducting electrons in the island, respectively. 

The second half of the parameter range e < Q g < 2e can be treated analogously. 
The tunneling now connects the states n = 2 and 1. The symmetry implies T eo ^ oe (Q g ) = 
Y eo l oe (2e—Q g ). This leads to the following picture for T < T cr : For — Q cr < Q g < Q cr (T) < e 
the system is in the even state, then it switches to the odd state where it stays for 
Qcr < Q g < 2e — Q cr , and periodic beyond. 

The switching point is temperature dependent. For A^ < Aj we find for T < Ti = 
2A t /\nN s 

Q cr {T) = ^+ e --^-\nN{T) (11) 

whereas for T} < T < T cr 

C p CT 

Qcr(T) = -(A, + A,) + | - — lniV(T) . (12) 

For A; > Aj the result ( pi]) with N(T) = N$ is valid at any temperature T < T cr . 

The expressions for Q cr were derived for Aj < Ec = e 2 /2C. In the opposite case, 
Aj > for T < Tq tunneling of 2 electrons becomes more favourable and the two-state 
model is not sufficient anymore. Here we defined 

J (A, + A,-£ C )/ In iV(T) for A* < A, - E c , , 

10 \ 2(A t -E c )/ In N s for A, > A, - E c [ ' 

On the other hand, for T > Tq the results given above remain valid. 

In the case where the lead electrode is normal A; = we reproduce the results of Refs. 
|T|-[3[], which were derived from considerations of the free energy of the superconducting island 
only. If the electrode is superconducting the energy of excitations in the lead e\ > A t needs 
to be taken into account, since the transition rates depend on the total energy difference. In 
contrast to the excitation which is trapped in the island, the excitations in the lead diffuse 
away on a time scale fast in comparison to inverse tunneling rates. Hence, their energy 
cannot be regained. This behavior involving different time scales is not accounted for in an 
equilibrium free energy description. 

The analysis presented above can also be extended to describe even-odd effects in NSN 
and SSS transistors. In this system the charging energy depends also on the transport voltage 
and on the number of electrons transported through the transistor. The total capacitance 
of the island C = C g + C\ + C r defines Ec = e 2 /2C. The energy differences for tunneling in 
the left and right junctions are 

IW^- '"'^' (14) 
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For definiteness we assume G\ = C r and eU tr < 2Ec, and we define Q tr = CU tr - It is again 
sufficient to consider only one even and one odd state. The transitions in this system are 
described by a master equation with the rates 

r eo -> Tf + T e r ° (15) 

which depend on the energy differences 5Ei and 5E r . They again are the sum of 5E c hj/ r 
and the energies of the excitations created in the island and electrodes. A similar relation 
holds for the transition from the odd state to the even state. After the substitution (p~5|) the 
stationary solutions for the occupation probabilities are given by eq. (H), and the cross-over 
gate voltage {Q g } cr = Qcr follows from (|). It depends on the transport voltage, which 
opens the possibility to tune the cross-over condition. The effect becomes visible only for 
T < eU tr - Then the two dominant tunneling processes are tunneling in the left junction 
with rate Tf followed by a tunneling process in the right junction T° e . Both balance each 
other in equilibrium. The cross-over temperature is defined by the condition Q cr (T) w e/2. 
For the SSS transistor (with A; = Aj) it is 

T cr = [2A-eU tr /2}/\nN s . (16) 

Next we evaluate the current due to single electron tunneling through the NSN or SSS 
transistors. For high temperatures T > T cr the current is / = e[F^T° e — r^T° e ]/T s , 
where the rates are given in (^|). It has maxima at the points Q g = e/2 + ne. At lower 
temperature T < T cr the current is I = e^ r . It is observable in the window Q cr {T) < Q g < 
| + AC/e + Qtr/2 and exponentially small outside. In NSN transistors it reduces to 

J=l/[2ei2 t JV < (0)y < ] (17) 

A second current peak or plateau exists in the window 3e/2 — AC/e — Qtr/2 < Q g < 2e — Q cr . 
Both plateaus create a double structure which repeats 2e-periodically. For A + eU tr /2 > Eq 
these plateaus merge forming a 2e-periodic single plateau structure. Note that the current 
(19) is much larger (by a factor 10 2 for parameters of ref. [6]) than the cotunneling current 
[4]. It can explain the presence of a constant current 80 / A of "unknown origin" detected in 
Ref. @. 

The occupation probabilities W e and W , which are regulated by the single electron 
tunneling processes, also influence the supercurrent through SSS transistors and the Andreev 
reflection in NSN systems. If the SSS transistors is in the even state Cooper pairs can tunnel 
through the system at Q g = ±e, ±3e, ... for small transport voltages. This leads to a set of 
2e-periodic sharp resonant peaks in the I-V curves. The amplitude of this peak is reduced 
below the T = result I e = I (T = 0)W e (e). Analogously, if the system is in the odd state 
Cooper pair tunneling through the system occurs at Q g = 0, ±2e, with current peak 
amplitude I Q = I(T = 0)W o (0). These probabilities are 



W ^ e ) e -2(A-E c )/T + y Ns 



e -(2A+E c )/T 
W °M = e -(2A-E c )/T + 1 / Ns 



(18) 
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At T > T* = 2 A/ hx(Ns) both current peaks are equal in height and the behaviour of the 
system is e-periodic. At lower temperature the two peaks are distinct. A behaviour of this 
type has recently been observed in experiments MM, the T-dependence is in qualitative 
agreement with the results of Ref @. 

If we consider NSN transistors we have to distinguish two cases. For A < Ec the 
single electron tunneling is crucial as discussed above. In the other limit A > Ec the 
mechanism of Andreev reflection transferring 2 electrons becomes important [IU|,n|. Close 
to Qg — ± e ; ±3e, ... at T = the shape of the current resonance 



Ires 

has been calculated in Ref 



G(5U g ,U tr )(u tr 



AC 2 

a 

U tr C 2 



(SU g 
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(19) 

for tunneling between even states. At finite temperature, 
due to single electron tunneling, there exists a finite probability for the odd state. This leads 
to an additional set of peaks due to 2e tunneling between two odd states. The amplitudes 
of the even-even and odd-odd current resonances then are I ee = I res W e {U g — 2Ec, U tr ) for 



Q g ~ e and I 



I res Wo(U g ,Ui 
W e (Q g p 



W Q (Q g « 0) 



tr, 



for Q g p; 



; 0, where 

)/ T + l/[iV (T) cosh 



CT 



-(A-E-)/T + 1/[N Q (T) COSh 



eSQg 
CT 



-{A+E c -e\U g \)/T 



e -(A-E c +e\U 9 \)/T + 1/[N (T) COSh 



elltr 
2T 



(20) 



(21) 



Here 5Q g = Q g 



E c ±eU tr /2. The two e-shifted peaks acquire equal height and the 



picture becomes e-periodic above a cross-over temperature TjJ = [A + E^]/ In N (T) > T cr . 
The presence of the odd peaks I 00 has been clearly demonstrated in recent experiments ||. 
For metals with a short elastic mean free path I at low temperatures and transport 



voltages the effective conductance G in eq. (|19f) near the maximum of the current is 



G(0,U tr 



12d 



A s 



e 2 lR 2 p 2 F S A 2 - (E, 



c, 



arctan 



A + E t 



c 



\ A -25, 



c 



(22) 



Here S is the junction area, d the size of the normal leads which have resistance small 
compared to R t . In contrast to the result quoted in Ref. fl3|], the conductance 
finite at T = eU tr = 0. For larger temperature or voltage 6 



remains 

max(T, eUtr/2) > E d = 
vpl/3d 2 the conductance G is smaller than (|22|) by a factor of order (E d /5y/ 2 and roughly 
agrees with the result of Ref. [0] in this limit. The result ( p^ ) is consistent with the height 
of the experimentally observed current peaks @,|] . 

In conclusion, we have developed a theory of parity effects in small superconducting 
islands by analyzing the rates of electron tunneling. Nonequilibrium excitations in the 
superconducting leads generated by the tunneling processes turn out to be essential. Our 
theory explains a large number of recent experiments in the electron box or transistors and 
can easily be extended to more general situations. 
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Figure Captions : 



Fig. 1 : Energy of different states in the NS electron box. Shown are the regular single 
electron tunneling transitions T ± and the transition 7 due to the tunneling of the odd 
electron. The state with two excitations and energy E > 2A relaxes quickly to the ground 
state by quasiparticle recombination. 
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